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We study charmonium correlators at finite temperature using quenched lattice 
QCD simulations. Two analysis procedures are applied to extract information 
on the spectral function: the maximum entropy method, and the x 2 fit analyses 
including the constrained curve fitting. We focus on the low energy structure of 
the spectral function by applying the smearing technique. We first discuss the 
applicability of these methods to the finite temperatures by analyzing the data 
at T = with restricted numbers of degrees of freedom. Then we apply these 
methods to the correlators at T > 0. We find no indication of mass shift and finite 
width for the charmonium states at T ~ 0.9T C . The results at T ~ 1.1T C imply 
that bound-state-like structures may survive even above T c . 



1. Introduction 

Investigation of the QCD phase transition at finite temperature and density 
is one of the most important subjects in hadron physics. To create the quark 
gluon plasma (QGP) phase, in which quarks and gluons are deconfined 
and the chiral symmetry is restored, the high energy heavy ion collision 
experiments have been performed and the RHIC experiment at BNL is 
in progress. The J/yj suppression has been considered as one of the most 
important signals of formation of QGP. However, theoretical understanding 
of charmonium at finite temperature and density is not sufficient to draw 
a conclusive scenario. 

We investigate the charmonium properties at finite temperature using 
lattice QCD simulations which enable us to incorporate the nonpertur- 
bative effects of QCD. Observations of spatial correlations between quark 
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and antiquark suggested nontrivial strong correlations persist above T c ' . 
These imply that bound-state-like structures may survive even above T c . 

In present work, we analyze the correlators by reconstructing the spec- 
tral function from lattice data using the maximum entropy method 3 and 
the x 2 fit analysis 4 ' 5 . The latter includes the constrained curve fitting 
6 in addition to the standard x 2 fit. To circumvent severe restriction of 
degrees of freedom in Euclidean temporal direction, we employ anisotropic 
lattices on which temporal lattice spacing is finer than the spatial one (for 
our setup of the anisotropic lattices, see Ref. 7 ). At this stage, the simu- 
lations are performed in the quenched approximation (without dynamical 
quark effects). To enhance the contribution from low energy part in the 
spectral function, we apply the smearing technique by spatially extending 
the meson operators. The details of this work has been presented in Ref. 5 . 

2. Analysis procedures 

The spectral function contains direct information on the structure of the 
mesonic correlator, C(t) = ^2-(0(x,t)0^ (0,0)). Recently the maximum 
entropy method (MEM) has been successfully applied to the extraction of 
the spectral function from lattice data at T — without assuming specific 
form 3 . However, its applicability to problems at T > is not straightfor- 
ward, because available physical range of C(t) as well as number of data 
points are severely limited. In such a circumstance, we need in some way 
to check the applicability of the method. Once a presumable form of the 
spectral function is known by MEM, other analysis procedures, such as the 
standard \ 2 fit may provide more quantitative information 4,5 . We there- 
fore propose to use the x 2 fit analysis in combination with MEM 5 . As 
a more sophisticated alternative to the standard x 2 fit, we also apply the 
constrained curve fitting 6 assuming several peaks with finite widths. In 
this method, one can add prior knowledge 1 for the spectral function and 
thus can treat functional forms with numbers of parameters (even) more 
than the number of data points. 

To verify the applicability of these methods at T > 0, we require that 
they correctly work for the correlators at T = with restricted numbers of 
data points corresponding to the cases at T > 0. As we will see in the next 
section, MEM applied to the correlators of the point operators does not 
meet this condition with present level of statistics (O(500) configurations). 

We therefore apply the smearing technique which enhances the low en- 
ergy part of the spectral function. For the correlators smeared with the 
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Figure 1. The spectral functions determined by MEM for the point correlator at T = 
in PS channel. 



wave function at T = 0, both the procedures work well in the low energy 
region. Although the features of the collective mode such as the mass and 
width are unchanged by the smearing, the smearing can produce a mimic 
peak structure in the spectral function 8 . This possibility is examined by 
observing the correlators smeared with a narrower function {half- smeared 
correlators). The dependence of extracted spectral function on the smear- 
ing function is compared with that of correlators composed of free quarks. 



3. Numerical results 

The numerical simulation is performed on lattices with the spatial lattice 
cutoff a^ 1 ~ 2 GeV and the anisotropy a a /a T = 4 in the quenched ap- 
proximation 7 . The lattice sizes are 20 3 x N t , where N t = 160 (T ~ 0), 32 
(T ~ 0.9T C ), and 26 (T ~ 1.1T C ). The numbers of configurations arc 500 at 
T ~ and 1000 at T > 0. N t — 28 roughly corresponds to the transition 
temperature. The quark field is described by the 0(a) improved Wilson 
action with the tree- level tadpole- improvement. The hopping parameter is 
chosen so that the charmonium spectrum is roughly reproduced. 



3.1. Results at T — 

Let us start with the analysis of the point correlators at T = 0. We apply 
MEM to them with restricted numbers of degrees of freedom. The results 
with two types of such restrictions are displayed in Fig. 1. The left panel 
shows the dependence of the result on t max , the maximum t of the correlator 
used in the analysis. This case corresponds to the situation at T > 0. MEM 
fails to reproduce even the lowest peak for t max < 16. The right panel shows 
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Figure 2. The spectral functions from the \ 2 fit analysis for the smeared correlators in 
the vector channel at T ~ 1.1T C . 

the results when one alternatively skips several time slices in the analysis. 
This case corresponds to the coarsening of the temporal lattice spacing. 
Even for t sep — 8 for which the number of data points is 6, MEM at 
least reproduces the correct lowest peak position while the resolution is not 
enough. These results indicate that the physical region of the correlator 
as well as the number of the degrees of freedom is important for MEM 
to work correctly. The required region of C(t) in the above analysis is 
tmax > O(0. 5fm), which is not fulfilled around T ~ T c . 

This is why we apply the smearing technique. The results of MEM for 
the smeared correlators are stable under the above two kinds of restriction 
for t m ax of interest; at least the lowest peak position is correctly reproduced 
while the resolution becomes worse as t max decreases. Since the smearing 
suppresses the high frequency region of the spectral function, we focus only 
on the lowest peak of the spectral function in the following analysis. 

The results of the \ 2 fit and the constrained curve fitting for the smeared 
correlators are consistent with the results of MEM and stable under the 
same restriction of t max . In the following, we concentrate on the results of 
the smeared correlators. 

3.2. Results at T > 

At T ~ 0.9T C , the result of MEM for the smeared correlator indicates that 
the low energy structure of the spectral function is well represented by a 
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strong peak corresponding to a meson state. We then apply the % 2 fit anal- 
ysis assuming double-pole, single relativistic Breit-Wigner (BW) type form 
4 , and BW+pole form (the pole function is to subtract the excited state 
contribution). The analysis indicates that the lowest peak of the spectral 
function is well represented by the pole form with the width consistent with 
zero. The mass is almost the same as at T = 0: no mass shift is observed 
at T ~ 0.9T C . The result of the constrained curve fitting is consistent with 
that of the standard x 2 fit. 

At T ~ 1.1 T c the result of MEM still exhibits a peak structure around 
the same energy region as at T < T c . Therefore we perform the same kind 
of the x 2 fit analysis as at T < T c . The fits to the 2-pole ansatz and to 
the Breit-Wigner type ansatze give inconsistent results, and the latter fits 
indicate that the spectral function has a peak with almost the same mass 
as at T < T c and the width of order of 200 MeV, as shown in Fig. 2. The 
result depends on the smearing function only slightly, in contrast to the 
free quark case, and hence we conclude that this is a physically significant 
structure. Similar result is obtained with the constrained curve fitting. 
These results are in accord with our earlier analysis of the spatial correlation 
between quark and antiquark 2 and imply that quasi-stable bound-state- 
like structures exist at T ~ 1.1T C . Similar results have also been reported 
by other groups 9 ' 10 . 

The simulation has been done on NEC SX-5 at Research Center for Nu- 
clear Physics, Osaka University and Hitachi SR8000 at KEK (High Energy 
Accelerator Research Organization). H. M. and T. U. are supported by 
Japan Society for the Promotion of Science for Young Scientists. 
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